Genome sequences of Tropheus moorii and Petrochromis trewavasae, two eco-morphologically divergent cichlid fishes endemic to Lake Tanganyika.

With more than 1000 species, East African cichlid fishes represent the fastest and most species-rich vertebrate radiation known, providing an ideal model to tackle molecular mechanisms underlying recurrent adaptive diversification. We add high-quality genome reconstructions for two phylogenetic key species of a lineage that diverged about ~ 3-9 million years ago (mya), representing the earliest split of the so-called modern haplochromines that seeded additional radiations such as those in Lake Malawi and Victoria. Along with the annotated genomes we analysed discriminating genomic features of the study species, each representing an extreme trophic morphology, one being an algae browser and the other an algae grazer. The genomes of Tropheus moorii (TM) and Petrochromis trewavasae (PT) comprise 911 and 918 Mbp with 40,300 and 39,600 predicted genes, respectively. Our DNA sequence data are based on 5 and 6 individuals of TM and PT, and the transcriptomic sequences of one individual per species and sex, respectively. Concerning variation, on average we observed 1 variant per 220 bp (interspecific), and 1 variant per 2540 bp (PT vs PT)/1561 bp (TM vs TM) (intraspecific). GO enrichment analysis of gene regions affected by variants revealed several candidates which may influence phenotype modifications related to facial and jaw morphology, such as genes belonging to the Hedgehog pathway (SHH, SMO, WNT9A) and the BMP and GLI families.

Asian horses deepen the MSY phylogeny.

Humans have shaped the population history of the horse ever since domestication about 5500 years ago. Comparative analyses of the Y chromosome can illuminate the paternal origin of modern horse breeds. This may also reveal different breeding strategies that led to the formation of extant breeds. Recently, a horse Y-chromosomal phylogeny of modern horses based on 1.46 Mb of the male-specific Y (MSY) was generated. We extended this dataset with 52 samples from five European, two American and seven Asian breeds. As in the previous study, almost all modern European horses fall into a crown group, connected via a few autochthonous Northern European lineages to the outgroup, the Przewalski's Horse. In total, we now distinguish 42 MSY haplotypes determined by 158 variants within domestic horses. Asian horses show much higher diversity than previously found in European breeds. The Asian breeds also introduce a deep split to the phylogeny, preliminarily dated to 5527 ± 872 years. We conclude that the deep splitting Asian Y haplotypes are remnants of a far more diverse ancient horse population, whose haplotypes were lost in other lineages.

Uncovering the genetic signature of quantitative trait evolution with replicated time series data.

The genetic architecture of adaptation in natural populations has not yet been resolved: it is not clear to what extent the spread of beneficial mutations (selective sweeps) or the response of many quantitative trait loci drive adaptation to environmental changes. Although much attention has been given to the genomic footprint of selective sweeps, the importance of selection on quantitative traits is still not well studied, as the associated genomic signature is extremely difficult to detect. We propose 'Evolve and Resequence' as a promising tool, to study polygenic adaptation of quantitative traits in evolving populations. Simulating replicated time series data we show that adaptation to a new intermediate trait optimum has three characteristic phases that are reflected on the genomic level: (1) directional frequency changes towards the new trait optimum, (2) plateauing of allele frequencies when the new trait optimum has been reached and (3) subsequent divergence between replicated trajectories ultimately leading to the loss or fixation of alleles while the trait value does not change. We explore these 3 phase characteristics for relevant population genetic parameters to provide expectations for various experimental evolution designs. Remarkably, over a broad range of parameters the trajectories of selected alleles display a pattern across replicates, which differs both from neutrality and directional selection. We conclude that replicated time series data from experimental evolution studies provide a promising framework to study polygenic adaptation from whole-genome population genetics data.

Spatially varying selection shapes life history clines among populations of Drosophila melanogaster from sub-Saharan Africa.

Clines in life history traits, presumably driven by spatially varying selection, are widespread. Major latitudinal clines have been observed, for example, in Drosophila melanogaster, an ancestrally tropical insect from Africa that has colonized temperate habitats on multiple continents. Yet, how geographic factors other than latitude, such as altitude or longitude, affect life history in this species remains poorly understood. Moreover, most previous work has been performed on derived European, American and Australian populations, but whether life history also varies predictably with geography in the ancestral Afro-tropical range has not been investigated systematically. Here, we have examined life history variation among populations of D. melanogaster from sub-Saharan Africa. Viability and reproductive diapause did not vary with geography, but body size increased with altitude, latitude and longitude. Early fecundity covaried positively with altitude and latitude, whereas lifespan showed the opposite trend. Examination of genetic variance-covariance matrices revealed geographic differentiation also in trade-off structure, and QST -FST analysis showed that life history differentiation among populations is likely shaped by selection. Together, our results suggest that geographic and/or climatic factors drive adaptive phenotypic differentiation among ancestral African populations and confirm the widely held notion that latitude and altitude represent parallel gradients.

Combining experimental evolution with next-generation sequencing: a powerful tool to study adaptation from standing genetic variation.

Evolve and resequence (E&R) is a new approach to investigate the genomic responses to selection during experimental evolution. By using whole genome sequencing of pools of individuals (Pool-Seq), this method can identify selected variants in controlled and replicable experimental settings. Reviewing the current state of the field, we show that E&R can be powerful enough to identify causative genes and possibly even single-nucleotide polymorphisms. We also discuss how the experimental design and the complexity of the trait could result in a large number of false positive candidates. We suggest experimental and analytical strategies to maximize the power of E&R to uncover the genotype-phenotype link and serve as an important research tool for a broad range of evolutionary questions.

A Drosophila laboratory evolution experiment points to low evolutionary potential under increased temperatures likely to be experienced in the future.

The ability to respond evolutionarily to increasing temperatures is important for survival of ectotherms in a changing climate. Recent studies suggest that upper thermal limits may be evolutionary constrained. We address this hypothesis in a laboratory evolution experiment, encompassing ecologically relevant thermal regimes. To examine the potential for species to respond to climate change, we exposed replicate populations of Drosophila melanogaster to increasing temperatures (0.3 °C every generation) for 20 generations, whereas corresponding replicate control populations were held at benign thermal conditions throughout the experiment. We hypothesized that replicate populations exposed to increasing temperatures would show increased resistance to warm and dry environments compared with replicate control populations. Contrasting replicate populations held at the two thermal regimes showed (i) an increase in desiccation resistance and a decline in heat knock-down resistance in replicate populations exposed to increasing temperatures, (ii) similar egg-to-adult viability and fecundity in replicate populations from the two thermal regimes, when assessed at high stressful temperatures and (iii) no difference in nucleotide diversity between thermal regimes. The limited scope for adaptive evolutionary responses shown in this study highlights the challenges faced by ectotherms under climate change.

Reproductive and post-reproductive life history of wild-caught Drosophila melanogaster under laboratory conditions.

The life history of the fruit fly (Drosophila melanogaster) is well understood, but fitness components are rarely measured by following single individuals over their lifetime, thereby limiting insights into lifetime reproductive success, reproductive senescence and post-reproductive lifespan. Moreover, most studies have examined long-established laboratory strains rather than freshly caught individuals and may thus be confounded by adaptation to laboratory culture, inbreeding or mutation accumulation. Here, we have followed the life histories of individual females from three recently caught, non-laboratory-adapted wild populations of D. melanogaster. Populations varied in a number of life-history traits, including ovariole number, fecundity, hatchability and lifespan. To describe individual patterns of age-specific fecundity, we developed a new model that allowed us to distinguish four phases during a female's life: a phase of reproductive maturation, followed by a period of linear and then exponential decline in fecundity and, finally, a post-ovipository period. Individual females exhibited clear-cut fecundity peaks, which contrasts with previous analyses, and post-peak levels of fecundity declined independently of how long females lived. Notably, females had a pronounced post-reproductive lifespan, which on average made up 40% of total lifespan. Post-reproductive lifespan did not differ among populations and was not correlated with reproductive fitness components, supporting the hypothesis that this period is a highly variable, random 'add-on' at the end of reproductive life rather than a correlate of selection on reproductive fitness. Most life-history traits were positively correlated, a pattern that might be due to genotype by environment interactions when wild flies are brought into a novel laboratory environment but that is unlikely explained by inbreeding or positive mutational covariance caused by mutation accumulation.

Multiple quantitative trait loci influence intra-specific variation in genital morphology between phylogenetically distinct lines of Drosophila montana.

The evolution of animal genitalia has gained renewed interest because of their potential roles during sexual selection and early stages of species formation. Although central to understanding the evolutionary process, knowledge of the genetic basis of natural variation in genital morphology is limited to a very few species. Using an outbred cross between phylogenetically distinct lines of Drosophila montana, we characterized quantitative trait loci (QTLs) affecting the size and shape of the distiphallus, a prominent part of the male intromittent organ. Our microsatellite-based linkage analysis shows that intra-specific variation in the distiphallus involves several QTLs of largely additive effect and that a highly significant QTL co-localizes with the same inversion where we have earlier localized a large QTL for a sexually selected courtship song trait. The latter indicates that inversions can play an important role in shaping the evolution of rapidly evolving traits with a potential influence on speciation.

Phylogeography, genetic structure and population divergence time of cheetahs in Africa and Asia: evidence for long-term geographic isolates.

The cheetah (Acinonyx jubatus) has been described as a species with low levels of genetic variation. This has been suggested to be the consequence of a demographic bottleneck 10 000-12 000 years ago (ya) and also led to the assumption that only small genetic differences exist between the described subspecies. However, analysing mitochondrial DNA and microsatellites in cheetah samples from most of the historic range of the species we found relatively deep phylogeographic breaks between some of the investigated populations, and most of the methods assessed divergence time estimates predating the postulated bottleneck. Mitochondrial DNA monophyly and overall levels of genetic differentiation support the distinctiveness of Northern-East African cheetahs (Acinonyx jubatus soemmeringii). Moreover, combining archaeozoological and contemporary samples, we show that Asiatic cheetahs (Acinonyx jubatus venaticus) are unambiguously separated from African subspecies. Divergence time estimates from mitochondrial and nuclear data place the split between Asiatic and Southern African cheetahs (Acinonyx jubatus jubatus) at 32 000-67 000 ya using an average mammalian microsatellite mutation rate and at 4700-44 000 ya employing human microsatellite mutation rates. Cheetahs are vulnerable to extinction globally and critically endangered in their Asiatic range, where the last 70-110 individuals survive only in Iran. We demonstrate that these extant Iranian cheetahs are an autochthonous monophyletic population and the last representatives of the Asiatic subspecies A. j. venaticus. We advocate that conservation strategies should consider the uncovered independent evolutionary histories of Asiatic and African cheetahs, as well as among some African subspecies. This would facilitate the dual conservation priorities of maintaining locally adapted ecotypes and genetic diversity.

A microsatellite linkage map for Drosophila montana shows large variation in recombination rates, and a courtship song trait maps to an area of low recombination.

Current advances in genetic analysis are opening up our knowledge of the genetics of species differences, but challenges remain, particularly for out-bred natural populations. We constructed a microsatellite-based linkage map for two out-bred lines of Drosophila montana derived from divergent populations by taking advantage of the Drosophila virilis genome and available cytological maps of both species. Although the placement of markers was quite consistent with cytological predictions, the map indicated large heterogeneity in recombination rates along chromosomes. We also performed a quantitative trait locus (QTL) analysis on a courtship song character (carrier frequency), which differs between populations and is subject to strong sexual selection. Linkage mapping yielded two significant QTLs, which explained 3% and 14% of the variation in carrier frequency, respectively. Interestingly, as in other recent studies of traits which can influence speciation, the strongest QTL mapped to a genomic region partly covered by an inversion polymorphism.

The impact of shared ancestral variation on hybrid male lethality--a 16 codon indel in the Drosophila simulans Lhr gene.

The Lethal hybrid rescue (Lhr) gene causes hybrid male lethality in crosses between Drosophila simulans and D. melanogaster. Lhr(2) is a D. simulans allele, which rescues hybrid males. It has been recently proposed that a 16 codon insertion, which distinguishes the D. melanogaster and the canonical D. simulans allele, and is lacking in Lhr(2), may be responsible for the functional divergence of D. melanogaster and D. simulans Lhr alleles. Here, we show that the Lhr(2) allele lacking the insertion represents an ancestral polymorphism segregating at a moderate frequency in D. simulans. Crosses of D. melanogaster females to males from two D. simulans strains carrying this deletion showed a severe deficiency of viable hybrid males. Our results suggest that the absence of this insertion alone is not sufficient to explain functional differences between D. melanogaster and D. simulans Lhr alleles.

Phylogeographic patterns in Drosophila montana.

The Drosophila virilis species group offers valuable opportunities for studying the roles of chromosomal re-arrangements and mating signals in speciation. The 13 species are divided into two subgroups, the montana and virilis 'phylads'. There is greater differentiation among species within the montana phylad in both karyotype and acoustic signals than exists among members of the virilis phylad. Drosophila montana is a divergent species which is included in the montana phylad. Here, we analyse the phylogeography of D. montana to provide a framework for understanding divergence of acoustic signals among populations. We analysed mitochondrial sequences corresponding to the cytochrome oxidase I and cytochrome oxidase II genes, as well as 16 microsatellite loci, from 108 lines of D. montana covering most of the species' range. The species shows a clear genetic differentiation between North American and Scandinavian populations. Microsatellite allele frequencies and mitochondrial DNA haplotypes gave significant FST values between populations from Canada, USA and Finland. A Bayesian analysis of population structure based on the microsatellite frequencies showed four genetically distinct groups, corresponding to these three populations plus a small sample from Japan. A network based on mitochondrial haplotypes showed two Finnish clades of very different shape and variability, and another clade with all sequences from North America and Japan. All D. montana populations showed evidence of demographic expansion but the patterns inferred by coalescent analysis differed between populations. The divergence times between Scandinavian and North American clades were estimated to range from 450,000 to 900,000 years with populations in Canada and the USA possibly representing descendants of different refugial populations. Long-term separation of D. montana populations could have provided the opportunity for differentiation observed in male signal traits, especially carrier frequency of the song, but relaxation of sexual selection during population expansion may have been necessary.

Silene tatarica microsatellites are frequently located in repetitive DNA.

The genomic distribution of microsatellites can be explained by DNA slippage, slippage like processes and base substitutions. Nevertheless, microsatellites are also frequently associated with repetitive DNA, raising the question of the relative contributions of these processes to microsatellite genesis. We show that in Silene tatarica about 50% of the microsatellites isolated by an enrichment cloning protocol are associated with repetitive DNA. Based on the flanking sequences, we distinguished seven different classes of repetitive DNA. PCR primers designed for the flanking sequences of an individual clone amplified a heterogeneous family of repetitive DNA. Despite considerable variation in the flanking sequence (pi = 0.108), the microsatellite repeats did not show any evidence for decay. Rather, we observed the emergence of a new repeat type that probably arose by mutation and was spread by replication slippage. In fact, a complete repeat type switch could be observed among the analysed clones. We propose that the analysis of microsatellite sequences embedded in repetitive DNA provides a hitherto largely unexplored tool to study microsatellite evolution.

Low abundance of Escherichia coli microsatellites is associated with an extremely low mutation rate.

It is widely assumed that microsatellites are generated by replication slippage, a mutation process specific to repetitive DNA. Consistent with their high mutation rate, microsatellites are highly abundant in most eukaryotic genomes. In Escherichia coli, however, microsatellites are rare and short despite the fact that a high microsatellite mutation rate was described. We show that this high microsatellite instability depends on the presence of the F-plasmid. E. coli cells lacking the F-plasmid have extremely low microsatellite mutation rates. This result provides a possible explanation for the genome-wide low density of microsatellites in E. coli. Furthermore, we show that the F-plasmid induced microsatellite instability is independent of the mismatch repair pathway.

Patterns of microsatellite variation through a transition zone of a chromosomal cline in Drosophila americana.

Chromosomal rearrangements have been considered as important barriers to gene flow and were often used in the delineation of species. The original taxonomic designation of Drosophila americana americana and Drosophila americana texana is based on the presence/absence of a centric fusion between the X- and fourth chromosomes. D. a. americana presents the derived fused state, whereas Drosophila a. texana presents the freely segregating ancestral state. The degree of genetic separation between the two chromosomal forms is still controversial, with different genetic markers yielding contrasting results even when the same populations were analyzed. Using 27 polymorphic microsatellites, we re-evaluated patterns of genetic differentiation between six D. americana populations sampled through a transition zone of both chromosomal forms in the central United States. Our results clearly reject a scenario of two differentiated species forming a hybrid zone in a region of parapatry and indicate that gene flow minimizes genome-wide differentiation associated with the two chromosomal arrangements.

Non-African origin of a local beneficial mutation in D. melanogaster.

It is well understood that the out-of-Africa habitat expansion of D. melanogaster was associated with the fixation of many beneficial mutations. Nevertheless, it is not clear yet whether these beneficial mutations segregated already in Africa or originated outside of Africa. In this article, we describe an ongoing selective sweep specific to one European population. One microsatellite allele has increased in a population from The Netherlands to a frequency of 18%, whereas it is virtually absent in 12 other European populations. The selective sweep resulted in a genomic region of more than 600 kb that is identical by descent. This is probably the first evidence of a beneficial mutation that has arisen outside of Africa and has resulted in a selective sweep localized in a population from The Netherlands.

A multilocus microsatellite phylogeny of the Drosophila virilis group.

We used a set of 48 polymorphic microsatellites derived from Drosophila virilis to infer phylogenetic relationships in the D. virilis clade. Consistent with previous studies, D. virilis and D. lummei were the most basal species of the group. Within the D. montana phylad, the phylogenetic relationship could not be resolved. Special attention was given to the differentiation between D. americana texana, D. americana americana and D. novamexicana. Significant differences between these three groups were detected by F(ST) analyses. Similarly, a model-based clustering method for multilocus genotype data also provided strong support for the presence of three differentiated groups. This genome-wide differentiation between D. americana texana and D. americana americana contrasts with previous analyses based on DNA sequence data.

World-wide survey of an Accord insertion and its association with DDT resistance in Drosophila melanogaster.

Previous work showed that insecticide resistance in Drosophila melanogaster is correlated with the insertion of an Accord-like element into the 5' region of the cytochrome P450 gene, Cyp6g1. Here, we study the distribution of the Accord-like element in 673 recently collected D. melanogaster lines from 34 world-wide populations. We also examine the extent of microsatellite variability along a 180-kilobase (kb) genomic region of chromosome II encompassing the resistance gene. We confirm a 100% correlation of the Accord insertion with insecticide resistance and a significant reduction in variability extending at least 20 kb downstream of the Cyp6g1 gene. The frequency of the Accord insertion differs significantly between East African (32-55%) and nonAfrican (85-100%) populations. This pattern is consistent with a selective sweep driving the Accord insertion close to fixation in nonAfrican populations as a result of the insecticide resistance phenotype it confers. This study confirms that hitchhiking mapping can be used to identify beneficial mutations in natural populations.